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STATUS OF THE EPIC/MOS CALIBRATION 
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Abstract 

The XMM-Newton observatory has the largest collect- 
ing area flown so far for an X-ray imaging system, result- 
ing in a very high sensitivity over a broad spectral range. 
In order to exploit fully these performances, a very accu- 
rate calibration of the XMM-Newton instruments is re- 
quired, and has led to an extensive ground and flight cali- 
bration program. We report here on the current status of 
the EPIC/MOS cameras calibrations, highlighting areas 
for which a reasonably good accuracy has been achieved, 
and noting points where further work is needed. 
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1. Introduction 

In order to achieve the highest collecting power ever de- 
ployed in X-ray astronomy, the XMM-Newton observa- 



tory carries three telescopes with identical mirrors ( Jansen 
et al. 2001). There is one CCD camera at the focus of 
each mirror, the ensemble of these three cameras mak- 
ing up the European Photon Imaging Camera instrument 



plane detectors. These hardware pieces, and the associated 
calibration items are : 

— the mirror, characterized by its on- axis effective area, 
vignetting, and point spread function, 

— the grating stack of the RGS (not present in the PN 
camera), characterized by its transmission function, 
and its azimuthal modulation, 

— the filter inside the EPIC camera, which can be either 
thin, medium or thick, characterized by its transmis- 
sion and its homogeneity, 

— the focal plane detectors, characterized by the quan- 
tum efficiency, charge transfer inefficiency, homogene- 
ity, redistribution matrix, energy scale of the CCDs, 
as well as their metrology. 

In addition to measuring the above characteristics, the 
calibration tasks also include the two problems of i) the 
background determination, essential for the study of ex- 
tended sources, and ii) pile-up treatment, essential for the 
study of strong point sources. 

Prior to launch, a large number of subsystem calibra- 
tions have been determined by ground measurements. The 
most outstanding are the mirrors characterization in UV 
and X-rays at the CSL and Panter facilities respectively, 
as well as the extensive CCDs measurements performed in 



the ORSAY/LURE synchrotron beams (e.g. Pigot et al 



(EPIC Turner et al. 2001). Two of the telescopes have 1999). If these ground results have given a solid basis for 



X-ray gratings mounted behind the mirrors, and for them 
the original flux is roughly equally shared between the 
EPIC foc al plane and the Reflect ion Grating Spectrome- 
ter (RGS, |den Herder et al. 2001 ); the two correspondin g 
EPIC cameras are made of MOS CCDs ({Short et al. 1998|) , 
and are the ones we will be discussing in this paper. The 
third telescope h as an unobstructed beam, and is equipped 
with PN CCDs ( fctriidcr et al. 200l|) . The cali bration sta- 
tus of this EPIC/PN camera is discussed by Briel et al 



(2002) in these proceedings. 



The overall response of the MOS cameras, which is the 
final product of interest for the observer when analysing 
his data, results from the combination of the calibration 



building the calibration data base, they had obviously first 
to be checked and refined with flight data. Furthermore, 
given the few percent accuracy goal set for the calibra- 
tions, a number of calibration aspects were not entirely, 
or not satisfactorily, covered by the ground measurements; 
this is the case for example of the CCD quantum efficiency 
at low energy. Some others, like the metrology or back- 
ground aspects, and obviously the monitoring of the CCD 
Charge Transfer Inefficiency degradation, can be studied 
with flight data only. This is the reason for the intensive 
flight calibration program which has followed the com- 
missioning phase, and for the routine calibration obser- 
vations which are scheduled (see list in Ehlc & Alticri 



of all tl ic individual hardware parts which arc encountered 
by the X-rays between the celestial source and the focal 



2001 ). A first re port on these in-fli ght calibration activi- 



ties was given by Lumb et al. (2000) . Although the exercise 
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is a rather drastic simplification of the complex calibration 
analysis, Table |l| gives a first cut estimate of the main con- 
tributor, ground or flight data, to each of the calibration 
items mentioned above. 



Table 1. First order calibration matrix, indicating the main 
source of data used for filling the Current Calibration 
Files. 





Item 


Ground 


Flight 


Telescope 


On axis area 


X 






Point Spread Function 




X 




Vignetting 




X 


Gratings 


Transmission 


X 






Azimuthal Modulation 


X 




Filters 


Transmission 


X 






Homogeneity 


X 




CCDs 


Quantum efficiency 


X 


X 




Charge Transf. Ineff. 




X 




Redistribution Matrix 


X 


X 




Homogeneity 


X 






Background 




X 




Metrology, astrometry 




X 




Timing 




X 



In this report, and after a brief section on the basic 
operating modes of the MOS cameras, we will mainly de- 
scribe the status of the items listed under the "flight data" 
column in Table [|. 

2. The EPIC/MOS cameras and modes 

The two EPIC/MOS cameras are made of 7 CCDs each, 
passively cooled to an operating temperature of —100° C . 
Each CCD is operated in a frame transfer mode, and has 
a 600 x 600 pixels imaging section, the pixel size being of 
40 fj,m, corresponding to l.l"on the sky. The CCDs are ar- 
ranged with a central CCD and 6 peripheral CCDs, as can 
be seen in Figure |l|. There are three basic imaging modes 
for the operation of the camera, and one timing mode. 
Images taken in these modes are displayed in Fig.jl]; the 
two top images correspond to MOS1 and MOS2 for the 
same field, and illustrate the fact the two MOS cameras 
are oriented at 90° from each other. For the timing mode, 
the short axis on the central chip image corresponds to 
a true spatial coordinate axis, while the long axis coordi- 
nate is here a measurement of time, and therefore carries 
no spatial information. Note that the peripheral CCDs 
are always operated in Full Frame, whatever the camera 
mode. 

3. Imaging characteristics 




Full Frame, MOS 1 & 2. Frame time : 2.6 s 




Large Window Small Window 

Frame time (sync) Frame time (sync) 

0.9 s central CCD 0.3 s central CCD 

2.7 s outer CCDs 2.7 s outer CCDs 




Timing 

Time resolution : 

1.5 ms central CCD 

Frame time 

2.6 s outer CCDs 



Figure 1. The four basic modes of the MOS cameras. 



3.1. Metrology and absolute astrometry 

Contrary to the PN camera which has a monolithic CCD, 
the seven CCDs of the MOS camera are mounted individ- 
ually, and it was necessary to measure their relative posi- 
tions, in X, Y, and rotation angle using celestial sources. 
The method for refining the MOS metrology involves iter- 
atively improving boresight and CCD offset values in the 
Current Calibration Files (CCFs) for both cameras, using 
known star fields. Starting with the most recent iteration 
of the CCFs calibrated event lists were produced using 
the MOS processing chain in the SAS, which were in turn 
used for generating cleaned, band limited, images suitable 
for source searching. The SAS source detection chain was 
then used to generate source lists. Using the maximum 
likelihood this source list was then used in conjunction 
with the USNO catalogue positions to i) remove gross off- 
sets in the central CCDs via shifts in the 3 boresight Euler 
angles, and ii) optimise the outer CCD displacements by 
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minimising the rms offset error between the detected po- 
sitions and catalogue positions for each CCD. 

This processing was performed on the Orion Molecular 
Cloud field, the NGC 2516 open cluster, and the Lockman 
Hole field. This has allowed the source detection error to 
be reduced to approx 0.4" and the rms offset errors be- 
tween detections and catalogue positions to approx 0.5"for 
MOS1 and 0.55"for MOS2. This accuracy, which is half 
the pixel size, is probably at the limit of this single field 
method and further improvements in the MOS metrology 
will have to come through statistical studies of large num- 
bers of detections over the lifetime of XMM-Newton. This 
work has now been fully reflected in the CCFs. 

It must be clear that these numbers characterize the 
internal focal plane metrology, and not the absolute point- 
ing accuracy reached in the absence of any positional ref- 
erence in the X-ray image. This last question has been 
addressed by Tedds fc Watson (2002) in these proceed- 
ings, who performed a statistical study on a very large 
number of fields; they have found that the absolute astro- 
metric error of the XMM-Newton coordinate frame has a 
distribution with a FWHM of ~ 3-4"per axis (R.A. and 
Dec), and 0.3°in roll angle. 

3.2. Point Spread Function 

The knowledge of the image shape (Point Spread Func- 
tion, PSF) of a point-like celestial source, and of the re- 
lated Encircled Energy Fraction (EEF) quantity, are obvi- 
ously of prime importance regarding the aspects of source 
detection and spectral extraction. The PSFs of the XMM- 
Newton mirrors have a rather complex shape, even for a 
source located on axis, as shown in Figure |^. Taking into 
account this complex shape is however unnecessary for 
most of the applications, in which the observer is simply 
interested in knowing the EEF for a given circular extrac- 
tion radius, or in comparing the observed radial profile to 
the expected one for a point-like source. 



For these purposes, Ghizzardi & Molendi (2002) have 



looked for a simple analytical representation of the PSF. 
They showed that a simple king profile, with only two 
parameters, represents the data very well. An example of 
fit is given in Figure 0. This study was done on a large 
set of data, up to 12 arcmin off-axis position. For each 
selected source, the core radius and the slope parameter 
were determined for seven energy intervals, when allowed 



by the statistics. Ghizzardi & Molendi (2002) found that 
these parameters have a simple linear dependence on the 
energy E and off-axis angle 9 of the incoming photons. 
They also provide a range of validity, both in E and 9, 
for this modeling, based on the availability of adequate 
point source observations. It must be noted in particular 
that the high energy and high off-axis angle regions are 
excluded from the range of application of this model. This 
parametrisation has been fully implemented in the CCFs. 





Figure 2. The PSF of MOS 1 (left) and MOS 2 (right) cam- 
eras. The images are 110" wide. 
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Figure 3. Example of a PSF fit, at 1 keV for a source being 
off- axis by 12* 
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Figure 4- Encircled Energy Fractions, for an on-axis 
source, in a ring with limiting radii of 5" and 60" . 



Although the precise core shapes look quite different 
for the two MOS cameras (Fig. ||), their encircled energy, 
at least on axis, is not drastically different. For example, 
Figure || shows the on-axis EEF for a ring with inner and 
outer radii of 5 and 60"respectively, typical of extraction 
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regions for piled-up sources. The difference in EEF for the 
two cameras, if not considered properly, would translate 
into a difference of about ~ 0.03 on the slope of a power 
law spectrum. 

3.3. Vignetting 



Compared to the PN camera, the MOS telescopes vi- 
gnetting is complicated by the presence of the gratings 
arrays. They introduced an azimuthal variation, as well 
as an additional energy dependence above about ~ 5 keV 
( Erd 2000 ). The in flight measurements were performed 
mainly with observations of G21.5-0.9, on axis and at 
10'and 12'off-axis positions. Although the variations with 
energy are not completely understood and are still un- 
der investigation, the uncertainty on the vignetting cor- 
rections is estimated to be better than ~ 5 % up to 10'off- 
axis, and ~ 10 % at 12'. 



EPIC-UOS Efficiency (Single Camera + Filter + Telescope) 
Open (Black); Thin niter (Red); Medium Filter (Green); Thick niter (Blue) 



4. Spectral response 

4.1. Effective area 

The EPIC/MOS effective displayed in Figure ||, 

combines the mirror area, the gratings and filter trans- 
missions, and the CCD quantum efficiency Starting from 
ground calibration data for all these quantities, Sembay 
et al. (2002) have succeeded in working out the corrections 
necessary to minimize the residuals obtained when fitting 
on-axis sources with simple spectra, such as BL Lacs. The 
main changes to the ground calibration data made in this 
work were on the Au edge for the mirror response, and 
the CCD quantum efficiency below one keV, a quantity 
which was poorly determined from the synchrotron mea- 
surements. Figure ^, from Sembay et al. (2002)| , illustrates 
the quality of the fit obtained on a spectrum of Mrk 421, 
with more than 2 millions counts per camera. The resid- 
uals are below 5% over all instrumental edges, and on the 
full range from 0.2 to 10 keV. 

Regarding now the overall spectrum, the quality of the 
fit does obviously not ensure the correctness of the pa- 
rameters. It does however allow detailed cross-calibrations 
between the EPIC instruments to be performed, as well 
as between EPIC and the RGS within XMM-Newton, or 
between EPIC and other X-ray instruments. This is dis- 



cussed separately in these proceedings by Griffiths et al 
(2002)|, [Haberl et al. (2002)| , |dcn Herder et al. (2002) 



and Snowden et al. (2002). There is generally good agree- 
ment between the EPIC instruments however small dis- 
crepancies still remain, namely differences in ur values 



(by ~ 1 x 10 



211 



cm 



on the low energy side, and differ- 




Energy (keV) 



Figure 5. Single EPIC/MOS camera efficiency, for the dif- 
ferent possible filters. The absorption edges due to the line- 
of-sight material have been indicated. 



Bl LooerUe Object. Mkn 421 
MPS! (Block) M0S2 (Red), Medium niter, ZEB Source Counts f Camera 




channel energy (keV) 



ences of ~ 0.1 for spectral indices on the high energy side. 
Besides that, future work will also be dedicated to the low- 



Figure 6. Fit and residuals of Mrk ^21 spectrum, for the 
two MOS cameras. 



est energies, below 200 eV, for which unexpected features 
have sometimes been observed. 



4.2. Charge Transfer Inefficiency 

Although they have been designed to be highly resilient 
to radiation, the MOS CCDs were expected to experi- 
ence a slow degradation of their parallel Charge Transfer 
Inefficiency due to their constant irradiation by high en- 
ergy protons. By reducing the magnitude of charge packets 
which originates furthest from the read-out node, CTI has 
the effect of shifting energies toward lower values, and to 
broaden lines. The monitoring of the CTI degradation is 
performed via the measurement of Al, Mn-Ka and -K/3 
lines generated by an internal calibration source, which is 
turned on via a special position of the filter wheel (closed 
w.r.t. the sky). These "calibration source" observations 
have been regularly scheduled at the beginning of each 
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orbit when XMM-Newton is still in the outer parts of the 
radiation belts in a too high radiation environment to al- 
low meaningful astrophysical observations. 
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Figure 8. Redistribution function at 600 eV (left) and fit 
to SNR 1E0102. 2-7219 (right). The red curve corresponds 
to the ground determination of the rmf, the black curve is 
the shape needed to obtain a good fit to the SNR data. 



rays which are absorbed close to the surface layer (Short 



TIME ELAPSED SINCE LAUNCH. 10-DEC- i999 (DAYSj 



et al. 1998 



Figure 7. Parallel transfer losses, for the Al and Mn-Ka 
lines since launch. The large solar flares times are indi- 
cated by the vertical lines. 



The detailed study of the CTI evolution is presented 
in these proceedings. The CTI has 



by Bcnnic et al. (2002 



been observed to increase smoothly with time, with jumps 
coincident with large solar flares, as shown in Figure ^. 
Algorithms have been designed to make a first order cor- 
rection to the CTI effect, and are included in the SAS 5.3 
version. Because of its statistical nature however, it is not 



possible to fully recover the CTI degradation. Bennie et al 
(2002)|have shown that there is still a 21 % broadening 
left on the FWHM of the Mn-Ka line measured on the 
full CCD chip, when the SAS 5.3 correction is used. 

Beyond improving the CTI correction in the form cur- 
rently implemented in SAS 5.3, a more drastic improve- 
ment under investigation is to adopt the concept of a devi- 
ation map as done for the Chandra ACIS. This would al- 
low to account for column to column individualities which 
have been evidenced, and reduce the amount of unrecover- 
able damage to ~ 5 eV/year for the Mn-Ka line FWHM, 
which had a pre-launch value of 125 eV. It remains that 
the CTI has an irreducible line broadening effect, which 
has now to be taken into account in the redistribution 
function. 

4.3. Redistribution function 

The ground calibration, in particular with monochromatic 
beams, has shown that the redistribution function (rmf) 
is rather complex. In particular, below 1 keV it has a sig- 
nificant low-energy shoulder which becomes increasingly 
important when the incident energy decreases. This is be- 
lieved to be due to an incomplete charge collection for X- 



Sembay et al. (2002) have been carefully modelling 
the rmf, using sources with strong emission lines, such 
as SNR 1E0102. 2-7219. As is shown in Figure § there 
has been a significant increase in the surface loss compo- 
nent, which required modifications to the rmf w.r.t. the 
pre-launch shape. The rmf are currently handled together 
with the effective area has stand-alone files. As mentioned 
above, the CTI degradation implies that a time dependent 
rmf has to be considered now; the way to implement that 
is under investigation. 

Finally, we note here that a single redistribution func- 
tion (as well as effective area) is found to be adequate 
for describing the response of the MOS for all imaging 
modes; it is not however adequate for the timing mode, 
which spectral response has still to be worked out. 



5. Background 

The background in EPIC has been triggering extensive 
discussions and work inside the EPIC calibration team 
(and outside. . . ), starting right after launch. A lot of work 
has been devoted to characterize this background and to 
look for ways of reducing it to a minimum. The XMM User 
Handbook gives a pretty complete account of the different 
components of the background found in EPIC, MOS and 
PN, which we will not be repeated here. We simply recall 
the main features of this background, and give the lat- 



est developments, in particular those from Lumb (2002a) 



Lumb (2002b) 



It is now recognized that the EPIC background is made 
of 4 components, being : 

1. A low energy electronic noise, 

2. A transient soft proton "flare" induced, 

3. A quiet time high energy protons induced, 

4. The diffuse cosmic X-ray background. 
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The first of this list is occuring at low energy, and is 
due to electronic interferences. It is selectively removed by 
the standard SAS processing. 

The second background is the one causing the strong 
"brightening" of the image inside the field of view, as in 
Figure [j]. It is induced by soft energy protons, below a 
few hundred keV, that directly reach the detectors via 
scattering through the mirrors. It happens in an unpre- 
dictable way, and has a highly variable light curve when it 
is pres ent. Moreover, it has been shown that its spectral 



shape ]s also highly variable. Despite numerous attempts, 
it was impossible to find any signature of these events 
that would allow them to be distinguished, and hence fil- 
tered out, from real X-rays. In practice, and except for 
very bright sources, periods contaminated by such a back- 
ground have to be excluded from the analysis. 

The characterization of the third background is the 
one on which most of the progress has been done. It was 
demonstrated to have only secular variations, as expected 
for an origin being from interactions of high energy cosmic- 
rays inside and around the detectors. Its spectrum is made 
of two components : a continuous component, and a line 
component. The continuous component has a hard spec- 
trum, displayed on Figure ||, probably generated by the 
combination of fluctuations in the energy losses suffered 
by the cosmic-rays when they go through the CCDs, and of 
Compton electrons due to interaction of the gamma-rays 
produced in the material around the CCDs. As expected 
because of its origin, there is no evidence for a spatial vari- 
ation of this component in the full focal plane. The line 
component is made up fluorescent lines, particularly Al 
and Si. Contrary to the continuous component, the distri- 
bution of these fluorescent lines, as well as weaker ones as 
Au for example, have a strong spatial dependency. This 
is evidenced in Figure [lO] which are maps made in the 
Al and Si lines respectively. They were generated from 
the sum of basically all the PV and CAL phases obser- 
vations, with only SNR observations excluded because of 
their strong line emissions. Because of potential residuals 
from the on-axis sources, the central part of the central 
CCD should not be considered significant here; however 
all the other structures are significant, and correspond to 
variations of about 15 % from one area to the other. This 
is readily explained by the localisation of the CCDs, also 
depicted in Figure |l^ : the outer CCDs detect more Al 
than the central one, because of their closer proximity to 
the Aluminium camera housing, while the Si excess seen 
in the central CCD is explained as Si escape lines from 
the back of the outer CCDs, which are located above the 
central one. The same explanation holds for the excess on 
the edges of CCDs 2, 4 and 6. 

In order to cope with this variable background, partic- 
ularly for extended sources, while awaiting for the develop- 
ment of an adequate SAS task to model it, |Lumb (2002a 



posure time, over 400 ks, sufficient for most observations, 
ii) an homogeneous collection of observations all at high 
galactic latitudes, with no bright sources, and all with the 
thin filter position. On each of the fields, a source detec- 
tion task was applied and led to about 10 objects per field 
to be excluded, corresponding typically to objects with a 
flux brighter than 1 — 2 x 10~ 14 ergs cm _2 s -1 . In the fi- 
nal set built by co-adding all observations, the excision of 
these point sources led to only a small 10 — 15 % level of 



local depression. As evidenced in Figures 3 and 4 of Lumt 



(2002b) , this set gives results very similar to the one noted 
above regarding the fluorescent lines spatial dependency. 

Finally, by comparing the spectrum derived from this 
data set with a model of the particle induced background 



spectrum, Lumb (2002b) has found a reasonable power 



law index of 1.4 for the ext ragalactic background. The 



same index is also found by de Luca fc Molcndi (2002) 



with a different data set and using a somewhat different 



MOSl red MOS2 black 




channel energy (keV) 



Figure 9. Spectrum of the internal background, in MOS1 
and MOS2 




has prepared a background template file and has made it 
available. This data set is characterised by i) a long ex- 



FigurelO. Spatial variations of the Al (left) and Si-Ka 
(right) lines of the internal background. 
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approach to model the particle background modelisation. 
The absolute normalisation however has still a rather large 
uncertainty. 

6. Areas in progress 

The above sections have stressed the areas for which most 
of the calibration work has been performed. In short, this 
is basically the understanding of the response to on-axis 
sources, observed in all imaging modes, the characteri- 
zation of the telescope vignetting, and of the background 
components. There are a number of calibration areas which 
have been less touched up to now, although this does not 
mean that they have particular problems. Some of them 
are briefly discussed below. 

6.1. Pile-up 

Although this is not strictly a calibration issue, it is im- 
portant to confirm the conditions for which pile-up has a 
sizeable effect on the spectra. Work is in progress on this 
area. Preliminary results confirm the theoretical analysis 



of Ballet 1999 which stresses that pile-up effects are min- 



imized by selection of isolated pixels (pattern 0). Await- 
ing for a potential correction of these pile-up effects, the 
observer has to exclude regions too heavily piled-up for 
spectral extractions. 

6.2. Off axis analysis 

At high energy, and apart from a slight difference in their 
depletion depth, the EPIC CCDs have homogeneous re- 
sponses. This is not the case below ~ 300 eV, for which 
ground measurements have shown evidence of large inho- 
mogeneities for some of the CCDs. While the high energy 
QE variations from one CCD to the other is implemented 
in the response functions, this is not the case yet of these 
inhomogeneities at very low energy within each CCD. One 
needs however to verify their presence in flight data, a task 
extremely difficult because of the interplay between redis- 
tribution and efficiency. Similarly, the thick filters have 
known inhomogeneities below the C edge which are not 
implemented yet. 

If a work as extensive as the one performed for on- 
axis sources has not been done yet, and will probably be 
impossible to do for the full focal plane given the obser- 
vation time it would request, once can say however that 
the analysis of "simple" extended sources, as clusters, has 
not revealed any significant problem for the overall off-axis 
response. 

6.3. Timing mode 
The timing functionalities of the timing mode ha ve been 



verified to be as expected. They are presented by Kuster 



et al. (2002) in these proceedings. The response function 



however is quite different from that of the imaging modes. 
This is due to the fact that rows are binned within the se- 
quencer program that reads out the CCDs. Also, contrary 
to imaging data, the timing data do not include any infor- 
mation on the energy in the pixels surrounding the X-ray 
event. This renders the tasks of noise cleaning and en- 
ergy determination extremely different from that of imag- 
ing data, and extremely difficult. This is still to be done 
properly, and at the time of writing, the spectra obtained 
in timing mode cannot be accurately modeled. 

7. Conclusions 

We have presented here a short summary of the status of 
the calibration of the MOS cameras of EPIC, as they are 
less than two years after the start of the initial calibra- 
tion phase. Most of the results have made, or are going 
to make, their way into version 5.3 of SAS. Apart from 
the references given here, the interested reader can find 
a lot of details on the calibration status in the technical 
notes maintained on the SOC web page. If the status for 
on-axis sources can be considered as quite good, we have 
also outlined a number of areas where work is still needed 
and will be done by the calibration team. 
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